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Light hadron masses with an 0(a 2 ) improved NNN action * 
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Meson and baryon masses in the light (u,d and s) sector are calculated using tadpole-improved gauge field and 
fermion actions. These are corrected to order 0(a 2 ) on the classical level using next-nearest-neighbour terms. 
The results, obtained at lattice spacings of 0.4 and 0.27fm, are compared to Wilson action calculations. 
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Improved actions hold great promise for mov- 
ing lattice field theory closer to phenomenology. 
Additional terms in the action that are nonlead- 
ing in powers of the lattice spacing a can be 
used to reduce discretization errors. These can 
be significantly suppressed further by incorporat- 
ing tadpole factors |l[ with the nonleading terms 
in the action. 

It is desirable to push improvement to the next- 
to-leading order in a. Alford et al proposed the 
so-called D234 action This action is cor- 

rected to 0(a 2 ) and tadpole improved. It is built 
by adding clover || and next-nearest-neighbour 
terms to the Wilson action. 

A simpler next-nearest-neighbour (NNN) ac- 
tion was considered some time ago Tree- 
level 0{a 2 ) improvement is achieved by employ- 
ing 2-link terms for the fermionic part and 6-link 
rectangles for the gauge field. We here use this 
NNN action with tadpole improvement to com- 
pute light hadron masses. Calculations with this 
action were also reported by Bock || and by 
Borigi and de Forcrand ||. 

The SU (3) gauge field part of the action is 
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The first term is the 4-link plaquette (pl) Wilson 
action and U r t are the planar 6-link rectangles 
(rt) which are sufficient to remove 0(a 2 ) errors at 
the classical level. The coefficient C r t = — 1/20U 2 
includes the tadpole factor 
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For the fermions NNN couplings in both the ki- 
netic and Wilson terms are used 10] 
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The Wilson action is recovered by replacing the 
coefficient 4/3 by 1 and dropping NNN terms. 

Calculations were carried out in quenched ap- 
proximation for the NNN and the Wilson action 
at matching physical lattice spacings a st deter- 
mined from the string tension |t]]. Two set of 
simulations were done at 0.4 and 0.27fm. The 
lattice parameters are shown in Tab. |]. 

Quark propagators were calculated for a range 
of k values using a stabilized biconjugate gradi- 
ent algorithm. Periodic boundary conditions were 
imposed on the quark fields in spatial directions 
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Table 1 



Lattice Details. Njj is the number of gauge configurations and a st is the lattice spacing determined from 
the string tension 0. k s is the hopping parameter corresponding to the strange quark mass. 



Action 


Lattice Njj 


(3 a st [fm] 




K 




NNN 


6 3 x 12 160 


6.25 0.4 


0.162, 0.165, 0.168, 0.171, 0.174 


0.166 


NNN 


8 3 x 14 60 


6.8 0.27 


0.148, 0.150, 0.152, 0.154, 0.156, 0.158 


0.1558 


Wilson 


6 3 x 12 160 


4.5 0.4 


0.189, 0.193, 0.197, 0.201, 0.205, 0.209, 0.213 


0.205 


Wilson 


8 3 x 14 90 


5.5 0.27 


0.164, 0.168, 0.172, 0.176, 0.180 


0.178 


Table 2 












Results of the calculations extrapolated to the limit M n — 0. The quantity J is defined in 


irj. 




Improved 




Wilson 


Exp. 




P = 6.25 


(3 = 6.8 


(3 = 4.5 


(3 = 5.5 




M p a p 


1.19(5) 


0.90(5) 


0.90(2) 


0.71(3) 




< 


648(27)MeV 


855(45)MeV 


858(15)MeV 


1085(46)MeV 




a st /a p 


1.31(5) 


1.17(7) 


1.73(3) 


1.50(6) 




J 


0.43(8) 


0.38(7) 


0.31(2) 


0.32(7) 




K/p 


0.65(2) 


0.65(4) 


0.61(1) 


0.61(4) 


0.64 


K*/p 


1.17(3) 


1.16(4) 


1.10(1) 


1.13(4) 


1.16 


Hp 


1.31(4) 


1.30(6) 


1.20(2) 


1.24(6) 


1.32 


N/p 


1.55(6) 


1.36(9) 


2.05(5) 


1.73(14) 


1.22 


A/N 


1.34(4) 


1.38(11) 


1.07(2) 


1.24(10) 


1.31 


t!/n 


1.15(2) 


1.20(4) 


1.05(1) 


1.10(5) 


1.27 


B/N 


1.23(3) 


1.32(4) 


1.09(1) 


1.15(7) 


1.40 


A/N 




1.17(4) 




1.08(5) 


1.19 


n~/N 


1.58(4) 


1.67(10) 


1.19(2) 


1.38(10) 


1.78 



but in the time direction a Dirichlet, or fixed, 
boundary condition was used. The source posi- 
tion was two time steps in from the boundary in 
all simulations. 

Meson and baryon correlators, constructed 
from standard local interpolating fields, were ob- 
tained for both local and smeared sinks with local 
sources. Gaussian smearing Bj was used. 

Masses were calculated using an analysis pro- 
cedure motivated by Bhattacharya et al || . This 
involves a combination of effective mass func- 
tions from local-local and local-smeared correla- 
tors. Except for the baryon spin-3/2 channel, it 
was found that compatible masses could be ob- 
tained using time averages starting 2 or 3 time 
steps away from the source. 

For each simulation the masses were extrapo- 
lated as a function of pion mass to the chiral limit 
M„ = 0. The choice of extrapolation function was 
either M = M + cM% or M = M + cM% + dM% 



depending on whether the coefficient d could be 
determined to be nonzero within the statistical 
errors. The errors in masses and in mass ratios 
were estimated using a bootstrap procedure with 
500 samples. 

In the strange quark sector we have fixed k s 
such that K*/K equals the experimentally ob- 
served 1.8. This was done by linear interpolation 
(or extrapolation) using the two k values from our 
set nearest to k s , see Tab. [[} 

Mass ratios extrapolated to the chiral limit are 
given in Tab. ||. The ratio M^/M p sets an overall 
scale of baryon masses relative to meson masses. 
This seems to be the quantity most effected by 
discretization errors. By 0.27fm the improved ac- 
tion results are fairly close to experiment and are 
compatible with Wilson action calculations done 
at small lattice spacing |ll| . 

An interesting way of looking at the relative 
quark mass dependence of pseudoscalar and vec- 
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Figure 1. The vector-pseudoscalar squared mass 
difference for the Wilson (solid symbols) and the 
NNN (open symbols) action for the (3 values of 
Tab. 0. 



Figure 2. Ratio of the lattice spacings determined 
by the string tension and the p-meson mass versus 
a p for the Wilson (solid symbols) and improved 
actions (open symbols). 



tor meson masses is to examine M v — Mp. Em- 
pirically, this quantity is approximately constant 
(s=s 0.5), a fact which has been linked to chiral 
symmetry pig ]. Therefore, its behaviour in lat- 
tice QCD may shed light on the role of chiral 
symmetry breaking Wilson-like terms in the ac- 
tion p3| . In Fig. [l]we plot the squared mass dif- 
ference in the limit of zero quark mass {M v = 0). 
At comparable lattice spacings the NNN results 



are much closer to M v 



Mp — const than the 



Wilson results. 

In Fig. U we have compiled some results for 
the ratio of the lattice spacing extracted from the 
string tension to the lattice spacing determined 
by the p-meson mass The improved 

action and Wilson action results show the same 
qualitative behaviour only shifted in lattice spac- 
ing by about a factor of 3. A remarkable feature 
seen in Fig. || is that different improved fermion 
actions exhibit a high degree of universality even 
in the non-scaling region. 

If simulations done with improved actions on 
coarse lattices are to be useful the results should 



extrapolate smoothly to the continuum limit. 
The ratio of nucleon to p-meson mass has been 
calculated a number of times. A sample of Wil- 
son action results [ 3||liyi7f| and results for tadpole 
improved actions 15 16[| are presented in Fig. |^. 
It is encouraging that our improved action values 
at 0.27 and 0.4fm are compatible with the Wilson 
action results at smaller lattice spacing. 

In this work we have seen that, compared to the 
Wilson action, the tadpole improved NNN action 
operates much closer to the continuum limit. Not 
suprisingly, the closeness of the mass spectrum to 
phenomenology, see Fig. 0, confirms this condu- 
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Figure 3. Ratio M^/M p versus a p for improved 
actions (open symbols) and the Wilson action 
(solid symbols). 



Figure 4. Level scheme of hadron masses from the 
NNN action at (3 = 6.8 and experimental values. 
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